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Statement  of  Problem 


The  overall  scientific  and  technical  effort  of  our  research  entails  research  in  scientific 
computations  concerning  the  control  of  the  flow  dynamics  of  key  semiconductor 
manufacturing  processes  that  operate  at  higher  flukes  of  nutrient  molecules  sustaining  film 
growth  than  accessible  with  conventional  processing  tools.  These  are  rapid  thermal 
organometallic  chemical  vapor  deposition  (RTOMCVD)  and  pulsed  Chemical  beam  epitaxy 
(PCBE).  The  modeling  is  closely  coordinated  with  hardware  design  and  process 
engineering  tasks  including  the  testing  of  novel  non-invasive  methods  of  real-time  process 
monitoring  and  control,  e.g.,  p-polarized  reflectance  spectroscopy.  The  development  of 
mathematical  models  has  been  initiated  that  accurately  describe  gas  flow  and  deposition 
process  for  RTOMCVD  reactors  in  various  operating  conditions  and  that  link  p-polarized 
reflectance  signal  to  the  composition  of  the  suface.  These  models  include  the  full 
compressible  gasdynamic  equations  which  incorporate  multiple  species  and  reactions 
between  the  species  coupled  with  boundary  conditions  that  describe  conditions  at  the  inlet, 
outlet  and  the  deposition  process.  Refinement  of  these  models  and  the  implementation  of 
closed-loop  process  control  are  objectives  of  continuing  research  funded  by  a  DoD  MURI 
on  Intelligent  materials  processing  and  electronics. 


Summary 


I.  Research  Accomplishments 


P-polarized  reflectance  spectroscopy  (PRS)  was  developed  as  a  method  of  real-time  pro¬ 
cess  monitoringC‘]-t^l.  The  relation  between  the  measured  reflectance  and  the  chemical 
kinetics  driving  the  growth  of  an  epitaxial  film  on  the  surface  of  a  substrate  has  been 
established  for  the  conditions  of  pulsed  chemical  beam  epitaxyt^l-t^f  The  precision  of  the 
monitoring  of  molecular  layer  growth  by  PRS  has  been  evaluatedt^f  Also,  limitations  to  the 
accuracy  of  PRS  monitoring  have  been  assessed.  The  design  of  the  RTOMCVD  reactor 
has  been  completed.  Experimental  validation  of  predictions  of  temperature  distributions  and 
flow  have  been  initiated.  We  have  applied  the  reduced  basis  method  for  control  of  the 
viscous  incompressible  gas  flow  dynamics  by  boundary  surface  controls  which  suggests 
that  real-time  simulation  and  control  can  be  done  by  lower-order  reduced  order  control 
systems. 


II.  Research  Accomplishments  and  Findings 


The  design  of  the  RTOMCVD  reactor  has  been  completed  and  numerical  simulations  of 
temperature  distributions  and  flow  in  this  reactor  have  been  performed.  We  have  chosen  a 


design  of  the  RTOMCVD  reactor  that  employs  a  channel  flow  across  the  substrate  wafer  - 
attached  to  a  rotating  susceptor  that  is  heated  by  radio-frequency  heating.  The  reactor  is 
comprised  of  (i)  an  inner  fused  silica  construction  that  shapes  the  flow  and  provides  for  a 
clean,  non-reactive  enclosure  for  chemical  vapor  deposition  process  and  (ii)  an  outer  steel 
shell  that  provides  for  containment  of  a  nitrogen  atmosphere.  This  outer  steel  shell  is  held  at 
the  chosen  pressure  with  controlled  pressure  equalization  between  the  inside  of  the  inner 
fused  silica  and  the  outer  steel  shells.  The  latter  is  presently  under  construction.  The  flow 
control  part  of  the  design  is  laid  out  for  four  process  gases  and  five  organometallic  (OlVl) 
source  vapors.  Provisions  are  made  for  pulsed  injection  of  the  process  gases  and  C^I 
source  vapors.  This  is  essential  for  the  operation  at  elevated  pressure  because  of  the  inverse 
relation  between  the  diffusion  constant  and  pressure,  requiring  OM  source  saturation  at  low 
pressure  combined  with  pulsed  injection  into  the  silica  flow  channel  by  plugs  of  high  pres¬ 
sure  carrier  gas.  Pulsed  operation  provides  also  for  instantaneous  probing  of  the  reaction 
kinetics  that  drives  the  heteroepitaxial  process  via  optical  real-time  process  monitoring  and 
mass  spectrometric  sampling  of  source  vapor  and  waste  vapor  compositions.  The  switch¬ 
ing  of  gases  is  designed  to  minimize  pressure  oscillations  in  the  reactor.  Applications  of 
RTOMCVD  will  focus  onto  heterostractures  employing  wide  bandgap  group  HI  phos¬ 
phides,  nitrides  and  phosphonitrides  in  the  context  of  optoelectronic  devicesf^l.  We  are 
assessing  at  present  by  infrared  imaging  the  temperature  distribution  for  a  specific 
prototype  fused  silica  shell  of  rectangular  cross  section  of  the  channel  for  atmospheric 
pressure  flow.  An  extension  of  this  assessments  to  a  scalable  multiwafer  reactor  is  planned. 


Initial  numerical  flow  simulations  on  such  a  RTOMCVD  chamber  have  been  performed  in 
order  to  investigate  the  effects  of  reactor  geometry  and  operating  conditions  on  flow  patters 
and  growth  rate  uniformity.  For  thermally  activated  process,  high  growth  rates  require 
operation  at  Iiigli  mass  flow  rate.  Because  this  reason  we  use  the  compressible  gas  dynaiu- 
ics  equations  to  model  the  flow  inside  a  2-D  axi-symmetric  reactor  chamber  approximating 
the  set  of  wafers  by  a  doughnut  shaped  ring.  We  carried  out  our  numerical  computations 
using  a  commercial  software  package  FIDAP,  for  the  finite  element  simulations  and  an  in- 
house  numerical  code  based  on  a  mixed  finite  element  method.  Our  numerical  simulations 
shows  that  under  the  current  design  of  the  horizontal  reactor  with  expanding  cross  section 
and  the  vertical  injection  of  the  inlet  gas  which  makes  the  ninety-degree  bend,  a  recircula¬ 
tion  region  exists  near  the  susceptor.  We  are  currently  studying  the  effect  of  the  operating 
conditions,  the  shape  of  the  inlet  region  and  thermal  distributions  on  reducing  this 
recirculation  cell.  These  and  other  simulation  results  provided  valuable  information  to  the 
design  of  an  optimal  reactor  chamber. 


Control  problems  that  involve  a  real-time  regulation  of  flow  and  growth  uniformity  in  the 
RTOMCVD  are  formidable  problems  to  solve.  The  state  equations  are  the  Navier-Stokes 
equations,  the.  geometry  is  often  complex  and  the  time  interval  involved  is  often  very  large. 
If  one  were  to  solve  such  problems  using  standard  finite  element  or  finite  difference  method 
the  resulting  system  is  prohibitively  large.  We  have  investigated  a  reduction  type  method 
based  on  the  so-called  reduced  basis  method,  which  over  comes  this  difficulty.  The 
reduced  basis  method  uses  basis  functions  which  are  closely  related  to  the  problem  that  is 
being  solved.  This  is  in  contrast  to  the  traditional  numerical  methods  such  as  finite 
difference  method  which  uses  grid  functions  as  basis  functions  or  finite  elements  method 
which  uses  piecewise  polynomials  for  this  purpose.  There  are  several  approaches  available 
for  the  selection  of  basis  functions.  We  tested  the  Lagrange  approach  that  uses  solutions  of 
the  problem  at  various  parameter  values  as  basis  functions.  The  major  advantage  of  this 
method  over  others  such  as  finite  element,  finite  difference  or  spectral  method  is  that  it  has 


fewer  degrees  of  basis  functions.  The  reduced  basis  which  we  developed  has  been  tested 
successfully  on  two  optimal  control  problems  in  steady  viscous  incompressible  flow  with 
boundary  control.  The  first  one  is  a  velocity  tracking  problem  in  cavity  flow  and  the  second 
one  is  vorticity  control  problems  in  channel  flow.  Currently  we  are  conducting  our  research 
on  developing  reduced-order  robust  feedback  methodologies  for  thermally  convected  flows 
in  channels  and  RTOMCVD  reactor. 


We  have  conducted  the  mathematical  studies  for  the  PRS  real-time  monitoring  technique. 
We  analyzed  the  PRS  signal  for  PCBE  epitaxial  film  growth  process  of  GaP/Si  hetero¬ 
structures.  We  have  carried  out  the  sensitivity  analysis  of  the  PRS  signal  with  respect  to  the 
modulation  of  the  substrate  temperature,  incident  angle  of  the  laser  beam  and  dielectric 
constants  of  the  film  and  substrate  layers,  and  measurement  noise.  The  relation  between  the 
measured  reflectance  and  the  activities  of  the  reactants  and  products  of  chemical  reactions 
that  feed  the  heteroepitaxial  growth  process  has  been  established  for  a  linearized 
approximation  of  the  complex  reflectivity.  Off-line  estimations  of  parameters;  the  incident 
angle  of  the  laser  beam,  complex  dielectric  constants  of  the  film  and  substrate  and  growth 
rate  per  pulsing  cycle  were  carried  out  successfully  using  the  least  square  method  for  data- 
fitting  problem,  under  homogeneous  growth  conditions.  We  are  currently  developing  on¬ 
line  recursive  algorithms  for  estimating  the  growth  rate  and  film  dielectric  constants  and 
feedback  control  laws  for  regulation  of  the  uniform  film  growth. 
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